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ABSTRACT 

Preliminary results are presented from a simple, single-antenna experiment designed to measure the 
all-sky radio spectrum between 100 and 200 MHz. The system used an internal comparison-switching 
scheme to reduce non-smooth instrumental contaminants in the measured spectrum to 75 mK. From 
the observations, we place an initial upper limit of 450 mK on the relative brightness temperature 
of the redshifted 21 cm contribution to the spectrum due to neutral hydrogen in the intergalactic 
medium (IGM) during the epoch of reionization, assuming a rapid transition to a fully ionized IGM 
at a redshift of 8. With refinement, this technique should be able to distinguish between slow and fast 
reionization scenarios. To constrain the duration of reionization to Az > 2, the systematic residuals 
in the measured spectrum must be reduced to 3 mK. 

Subject headings: Cosmology: Early Universe, Galaxies: Intergalactic Medium, Radio Lines: General 
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1. INTRODUCTION 

The transition period at the end of the cosmic "Dark 
Ages" is known as the epoch of reionization (FOR) . Dur- 
ing this epoch, radiation from the very first luminous 
sources — early stars, galaxies, and quasars — succeeded in 
ionizing the neutral hydrogen gas that had filled the in- 
tergalactic medium (IGM) since the recombination event 
following the Big Bang. Reionization marks a signifi- 
cant shift in the evolution of the Universe. For the first 
time, gravitationally-coUapsed objects exerted substan- 
tial feedback on their environments through electromag- 
netic radiation, initiating processes that have dominated 
the evolution of the visible baryonic Universe ever since. 
The epoch of reionization, therefore, can be considered a 
dividing line when the relatively simple evolution of the 
early Universe gave way to more complicated and more 
interconnected processes. Although the Dark Ages are 
known to end when the first luminous sources ionized 
the neutral hydrogen in the IGM, precisely when this 
transition occurred remains uncertain. 

The best existing constraints on the timing of the 
reionization epoch come from two sources: the cosmic 
microwave background (CMB) anisotropy and absorp- 
tion features in the spectra of high-redshift quasars. The 
amplitude of the observed temperature anisotropy in the 
CMB is affected by Thomson scattering due to elec- 
trons along the line of sight between the surface of last 
scattering and the detector, and thus, it is sensitive to 
the ionization history of the IGM through the electron 
column density. In addition, if there is sufficient opti- 
cal depth to CMB photons due to free electrons in the 
IGM after reionization, some of the angular anisotropy 
in the unpolarized intensity can be converted to polar- 
ized anisotropy. This produces a peak in the polariza- 
tion power spectrum at the angular scale size equivalent 
to the horizon at reionization with an amplitude pro - 
portional to the optical depth (jZaldarriaga et al.]|1997D . 



Measurements by the WMAP satellite of these effects 
indicate that the re dshift of reionization \s Zr ~ 11 ± 4 
(jSpergel et al.l[200l . assuming an instantaneous transi- 
tion. 

Lyman-a absorption by neutral hydrogen is visible in 
the spectra of many high-redshift quasars and, thus, of- 
fers the second currently feasible probe of the ionization 
history of the IGM. Continuum emission from quasars is 
redshifted as it travels through the expanding Universe 
to the observer. Neutral hydrogen along the line of sight 
creates absorption features in the continuum at wave- 
lengths corresponding to the local rest-frame wavelength 
of the Lyman-a line. Whereas CMB measurements place 
an integrated constraint on reionization, quasar absorp- 
tion line studies are capable of probing the ionization 
history in detail along the sight-lines. There is a signifi- 
cant limitation to this approach, however. The Lyman- 
a absorption saturates at very low fractions of neutral 
hydrogen (of order xhi ~ 10~^). Nevertheless, results 
from these studies have been quite successful and show 
that, while the IGM is highly ionized below z <Q (with 
typical Xhi ^5 10~^), a significant amount of neutral hy- 
drogen is present above, although precisel y how much re- 
mains unc lear fo iorgovski et al. 2001; Be cker et al.ll2001l: 
Fan et all 12003. 120031 : iWvithe fc Loebi 120041: iFan enil 
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The existing CMB and quasar absorption measure- 
ments are somewhat contradictory. Prior to these stud- 
ies, the reionization epoch was assumed generally to be 
quite brief, with the transition from an IGM filled with 
fully neutral hydrogen to an IGM filled with highly ion- 
ized hydrogen occurring very rapidly. These results, how- 
ever, open the possibility that the ionization history of 
the IG M may be more complicated than previously be- 
lieved (Haiman fc Holdej [2003l:[CCTl[200llSokasian et aH 
|2003; Madau et al. 200f . 

Direct observations of the 21 cm (1420 MHz) hy- 
perfine transition line of neutral hydrogen in the IGM 
during the reionization epoch would resolve the exist- 
ing uncertainties and reveal the evolving properties of 
the IGM. The redshifted 21 cm signal should appear as 
a faint, diffuse background in radio frequencies below 
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V < 200 MHz for redshifts above z > 6 (according to 
ly = 1420/ [1 + z] MHz). For diffuse gas in the high- 
redshift {z « 10) IGM, the expected unpolarized differ- 
ential brightness temperature of the redshifted 21 cm hue 
relative to the pervasive CMB is readily calculable from 
basic principles and is given by (.Zaldarriaga et al.ll2004l . 
their § 2) 



Hoean k ReeslHOTilScott fc Reeslll990tllliev et al ll2002L 



( i\ \( _ai 
^ \^ 0.02 ) \\n„^ 



Ts 



1/2 



mK, 



where S{9,z) is the local matter over-density, xhi{0,z) 
is the neutral fraction of hydrogen in the IGM, T^{z) = 
2.73 (1 + 2:) K is the temperature of CMB at the redshift 

of interest, Ts{9, z) is the spin temperature that describes 
the relative population of the ground and excited states 
of the hyperfine transition, and is the baryon density 
relative to the critical density, r2„i is the total matter 
density, and h specifies the Hubble constant according 
to Hq = 100 h km s~^ Mpc~^. From Equation [U we see 
that perturbations in the local density, spin temperature, 
and neutral fraction of hydrogen in the IGM would all be 
revealed as fluctuations in the brightness temperature of 
the observed redshifted 21 cm line. 

The differential brightness temperature of the red- 
shifted 21 cm line is very sensitive to the Hi spin tem- 
perature. When the spin temperature is greater than 
the CMB temperature, the line is visible in emission. 
For Ts ^ T^, the magnitude of the emission saturates to 
a maximum (redshift-dependent) brightness temperature 
that is about 25 to 35 mK for a mean-density, fully neu- 
tral IGM between redshifts 6 and 15, assuming a ACDM 
cosmology with 17,„ = 0.3, flA — 0.7, f2b — 0.04, and 
h = 0.7. At the other extreme, when the spin tempera- 
ture is very small and Ts <^T^, the line is visible in ab- 
sorption against the CMB with a potentially very large 
(and negative) relative brightness temperature. 

A number of factors are involved in predicting the typ- 
ical differential brightness temperature of the redshifted 
21 cm line as a function of redshift. In particular, the 
spin temperature must be treated in detail, including 
collisional coupling between the spin and kinetic tem- 
peratures of the gas, absorption of CMB photons, and 
heating by ultra-violet radiation from the first luminous 
sources. We direct the reader to IFurlanetto et al.l (|2006f ) 
for a good introduction to the topic. The results of 
several efforts to predict the evolution of the differen- 
tial brightness temperature of the redshifted 21 cm line 
have yielded predictions that are generally consistent 
in overall behavior , but vary highly in specific details 
(iMadau et al.lll9 97l:ISh "aver et aLlflOMlGnedin fc Shaved 
120041: IFurlanett o 2006). These models tend to agree 
that, for a finite period at sufficiently high redshifts 
(z > 20), the Hi hyperfine line should be seen in ab- 
sorption against the CMB, with relative brightness tem- 
peratures of up to I^Tfil < 100 mK. This is because 
the IGM initially cool s more rapidly than the CMB fol- 
lowing recombination (jShapiro et al.|[T99^ iMadau et al.l 
[1997). During this period, fiuctuations in the differ- 
ential brightness temperature of the redshifted 21 cm 
background should track t he underlying baryonic mat- 
ter density perturbations (jSunvaev fc ZeldovichI 119721 : 



20031: iLoeb & Zaldarriagal[200l iBarkana fc Loebl[2005h . 

Eventually, however, the models indicate that the ra- 
diation from the first generations of luminous sources 
will elevate the spin temperature of neutral hydrogen in 
the IGM above the CMB temperature and the redshifted 
21 cm line should be detected in emission with relative 
brightness temperatures up to the expected maximum 
values (of order 25 mK). Finally, during the reioniza- 
tion epoch, the neutral hydrogen becomes ionized, leav- 
ing little or no gas to produce the Hi emission, and 
the apparent differential brightness temperature of the 
redshifted 21 cm line falls to zero as reionization pro- 
gresses. As the gas is ionized, a unique pattern should 
be imprinted in the redshifted 21 cm signal that re- 
flects the processes responsible for the ionizing photons 
and that evolves wit h redshift as reion i zation progresses 
(IMadau et al.M l997t Pfozzi et al |[2000l: ICiardi fc Madaul 
[2003; Zaldarri agTet al.l 120041 : IFurlanetto et al.l 120041 ). 
The details of the specific timing, duration, and mag- 
nitude of these features remains highly variable between 
theoretical models due largely to uncertainties about the 
properties of the first luminous sources. 

Measuring the brightness temperature of the redshifted 
21 cm background could yield information about both 
the global and the local properties of the IGM. Deter- 
mining the average brightness temperature over a large 
solid angle as a function of redshift would eliminate 
any dependence on local density and temperature per- 
turbations and constrain the evolution of the product 
xhi{^ — T-f/Ts), where we use the bar to denote a spa- 
tial average. During the reionization epoch, it is, in gen- 
eral, believed to be a good approximation to assume that 
Ts ^ T^ and, therefore, that the brightness temperature 
is proportional directly to xri- Global constraints on 
the brightness temperature of the redshifted 21 cm line 
during the EOR, therefore, would directly constrain the 
neutral fraction of hydrogen in the IGM. Such constraints 
would provide a basic foundation for understanding the 
astrophysics of reionization by setting bounds on the du- 
ration of the epoch, as well as identifying unique features 
in the ionization history (for example if reionization oc- 
curred in two phases or all at once). They would also 
yield improvements in estimates of the optical depth to 
CMB photons and, thus, would help to break existing 
degeneracies in CMB measurements between the optical 
depth and prope rties of the primordia l matter density 
power spectrum (jTegmark et al.l [20061) . For these rea- 
sons, several efforts are underway to make precise mea- 
surements of the radio spectrum below v < 200 MHz 
(z > 6). In this paper, we report on the initial results 
of the Experiment to Detect the Global EOR Signature 
(EDGES). In §[21 we describe the specific approach used 
for EDGES to address the issue of separating the red- 
shifted 21 signal from the foreground emission. We then 
give an overview of the EDGES system in § |3l followed 
by the results of the first observing campaign with the 
system in § (H along with a discussion of the implications 
for future single- antenna measurements. 

2. METHOD 

In principle, the global brightness temperature mea- 
surement is much less complicated to perform than 
the detection of local perturbations in the redshifted 
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Fig. 1. — EDGES deployed at Mileura Station in Western Aus- 
tralia. The left panel shows the full antenna and ground screen in 
the foreground and the analog-to-digital conversion and data ac- 
quisition module in the background. The right panel is a close-up 
view of the amplifier and switching module connected directly to 
the antenna (through the balun). 

21 cm background (which will be attempted in the 
near future by the Mileura Widefield Array [MWA], 
Low Frequency Array [LOFAR] , Giant Metrewave Radio 
Telescope [GMRT], and Twenty-one Centimeter Array 
[21CMA]). Since the desired signal for the global mea- 
surement is the mean brightness temperature due to red- 
shifted 21 cm emission (or absorption) over the entire sky, 
there is no need for high angular resolution or imaging. A 
single antenna tuned to the appropriate frequencies could 
reach the required sensitivity (~ 25 mK) within only one 
hour of integration time, assuming a reasonable spectral 
resolution of ~ 1 MHz (equivalent to Az « 0.1 at z « 8). 
There is a fundamental complication with such an ex- 
periment, however, arising from the global nature of the 
signal. Since the expected redshifted 21 cm emission fills 
the entire sky, there is no ability to perform comparison 
switching between the target field and a blank field. The 
problem this causes is two-fold. First, it is difficult to 
separate the contribution to the measured spectrum due 
to the redshifted Hi signal from that of any other all-sky 
emission, including Galactic synchrotron and free- free ra- 
diation, the integrated contribution of cxtragalactic con- 
tinuum sources, or the CMB. Second, for similar reasons, 
it is difficult to avoid confusing any systematic effects in 
the measured spectrum due to the instrument or envi- 
ronment with received signal from the sky. The severity 
of these problems is exacerbated in single-antenna mea- 
surements by the intensity of the Galactic synchrotron 
emission. Unlike interferometric observations, a single 
antenna is sensitive to the large-scale emission from the 
Galaxy, providing a 200 to 10,000 K foreground in the 
measured spectrum. 

Determining the ~ 25 mK redshifted 21 cm contri- 
bution to the radio spectrum requires separating the 
signal from the foreground spectrum at better than 1 
part in 10,000. This can be accomplished by taking ad- 
vantage of the differences between the spectra of the 
Galactic and extragalactic foregrounds and the antici- 
pated reds hifted 21 cm contribution. As discussed by 
iShaver eta l. (1999), Galactic synchrotron emission is the 
dominate component of the astrophysical foregrounds be- 
low V < 200 MHz, accounting for all but approximatel y 
30 to 70 K of the foregrounds at 178 MHz ()Bridldll967t ). 



Its spectrum is very nearly a power-law given, in tem- 
perature units, by Tgai{v) ~ ^ where /3 « 2.5 is the 
spectral index. The spectral index is generally constant 
over the frequencies of interest (50 '^v< 200 MHz), al- 
though it is known to flatten with decreasing frequency 
due to self-absorption. The intensity of the synchrotron 
emission and the exact value of the spectral index depend 
on Galactic coordinate. The amplitude varies over an or- 
der of magnitude, between about 200 < Tgai < 10, 000 
K at 150 MHz (peaking toward the Galactic center and 
along the Galactic plane), while the spectral index has 
small variations of order ap « 0.1 dependent largely 
on Galactic latitude, with the steepest regions occur- 
ring at high Galactic latitudes. Free-free emission in the 
Galaxy and discrete Galactic and extragalactic contin- 
uum sources also have spectra that can be reasonably 
described by power-laws. The integrated flux from extra- 
galactic continuum sources is generally isotropic on large 
scales and accounts for the majority of the remaining 
power in the low- frequency radio spectrum, with free-free 
emission making up only about 1% of the total power. 
The combined spectrum due to the astrophysical fore- 
grounds is smooth and remains similar to a power-law 
profile. 

On the other hand, as the apparent differential bright- 
ness temperature of the redshifted 21 cm background 
transitions from T21 = mK at very high redshift to 
T21 « -100 mK before heating of the IGM by the first 
luminous sources, and then climbs to T21 ~ 25 mK at the 
beginning of the reionization epoch before falling back to 
T21 ~ mK at the end of reionization, the global mean 
redshifted 21 cm spectrum may contain up to three rel- 
atively sharp features between 50 < 1/ < 200 MHz that 
would not be represented well by a power-law profile. 
For the large solid angles of a single antenna beam, the 
mean redshifted 21 cm signal should vary little from one 
location to another on the sky. IFurlanettol ({2006) and 
iGnedin fc Shaverl ()2004[ ) have calculated example global 
mean redshifted 21 cm spectra for various assumptions 
of stellar formation histories. 

The specific approach employed with EDGES to ex- 
ploit these expected differences in spectral characteristics 
in order to overcome the difficulty in separating the fore- 
ground and signal contributions in the measured spec- 
trum is to limit the scope of the experiment to test for 
discontinuous features in the spectrum, since these would 
necessarily be due to the rapid transitions in the red- 
shifted 21 cm brightness temperature and not the spec- 
trally smooth foregrounds. In particular, the frequency 
response of the system is designed to test for fast reion- 
ization only (and not the transitions that might arise at 
higher redshifts from cooling and heating of the IGM). In 
the extreme case that the transition from a fully neutral 
to a fully ionized IGM was virtually instantaneous, such 
that XHi{zr) — 00, where Zr is the redshift of reioniza- 
tion, the contribution to the global spectrum at the fre- 
quencies corresponding to the reionization epoch would 
approach a step function. A sharp feature resembling a 
step function that is superimposed on the smooth power- 
law-like foreground spectrum should be relatively easy 
to identify. And if reionization were to progress more 
slowly, producing a smooth transition that spanned a 
large range of redshifts and many tens of MHz, a simple 
model could be used to set limits on the maximum rate 
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Fig. 2. — Example of redshifted 21 cm c ontribution (solid) to 
T^ky based on the model described in ^ 14.11 with AT21 = 25 mK, 
Zr = 8, and Az = 0.6. The residuals are also shown for a seventh- 
order polynomial fit to a simulated spectrum between 130 and 
190 MHz with (dash) and without (dot) the redshifted 21 cm con- 
tribution. The foreground contribution was modeled for the plot 
using 13 = 2.5, and Tgai(150 MHz) = 250 K. 

of the transition. 

In principle, a variety of such models could be devised 
to use in tests for the presence of a step feature in the 
radio spectrum due to a rapid reionization. A simple 
low-order polynomial fit to the measured spectrum would 
reveal such a discontinuous feature in the residual spec- 
trum after subtracting the fit and, thus, would be able 
to determine the redshift range of a rapid reionization. 
Figure [2] illustrates this approach by plotting a model 
(described in Section |4?T|) of the redshifted 21 cm contri- 
bution to the measured spectrum along with the residu- 
als after a seventh-order polynomial fit is removed from 
a simulated sky spectrum. This is the method used for 
the preliminary EDGES measurements. An advantage of 
this approach for global reionization experiments is that, 
given sufficient sensitivity, even a null result could still 
constrain xhi{z) and, thereby, distinguish between slow 
and fast reionization scenarios. 

3. EXPERIMENT DESIGN 

By focusing (at least initially) on confirming or rul- 
ing out a fast reionization scenario, the design of the 
EDGES system is able to be relatively simple. The pri- 
mary need is to reduce any instrumental or systematic 
contributions to the measured power spectrum that vary 
rapidly with frequency, since these could be confused 
with a sharp feature in the spectrum due to a fast reion- 
ization of the IGM. Such contributions could be due to 
terrestrial transmitters, reflections of receiver noise or 
sky noise from nearby objects, undesirable resonances 
within the electronics or the radio-frequency interference 
(RFI) enclosures, or spurious signals introduced by the 
digital sampling system. In this section, we provide an 
overview of the experimental design and setup, highlight- 
ing aspects that are relevant to reducing the effects of 
both the external and internal sources of systematic er- 
rors. Additional details on the analysis of systematic 
contributions and the hardware design can be found in 
the EDGES memorandum series^. 

^ http:/ /www. haystack. mit.cdu/ast/arrays/Edges/ 



3.1. Site Selection 

Some of the contributions to the systematic uncer- 
tainty listed above can be addressed by careful selec- 
tion of the observing site. Avoiding terrestrial transmit- 
ters (primarily from FM radio and television stations) 
is the most serious problem. Even at distances of hun- 
dreds or thousands of kilometers, tropospheric ducting 
and scattering (troposcatter), sporadic E propagation in 
the ionosphere, and reflections from meteors are all ca- 
pable of transferring a significant amount of power from 
Earth-based transmitters. The background produced by 
the integrated effect of many distant transmitters may 
have significant spectral structure above the expected 
redshifted 21 cm level. For example, a single, 100 kW 
FM radio station at 300 km from the observing site could 
produce up to a 100 K effective temperature in a 1 MHz 
channel due to troposcatter, or 100 mK due to meteor re- 
flections. Fortunately, these mechanisms of atmospheric 
propagation exhibit diurnal or transient behavior (as is 
the case for sporadic E propagation, tropospheric duct- 
ing, and meteor reflections) or require specific geometries 
for peak efficiency (as is the case for troposcatter) , mak- 
ing sensitive measurements possible at remote sites at 
least some of the time. 

Another concern is that local objects in the environ- 
ment will scatter both external noise and receiver noise, 
which will then be picked up by the system and corre- 
late with the original noise, causing sinusoidal ripples in 
the measured spectrum. We have estimated the magni- 
tude of the reflections of the Galactic foreground from 
objects like trees and mountains on the horizon where 
the antenna gain is reduced by a factor of 20 dB or 
more. As long as objects subtend solid angles under 
about 100 deg^, the spectrum will only be affected by 
a few parts per million (ppm). We have also considered 
the magnitude of noise originating from the receiver that 
will be returned by a nearby scatterer. Even if we assume 
that this noise is perfectly correlated with the internal re- 
ceiver noise, it will only produce ripples in the spectrum 
at the level of a few ppm provided that the object, like a 
tree subtending a few deg^, is more than ~ 100 m away 
or a larger object, subtending ^100 deg^ is more than 
~ 1 km away. 

Reflections of signals from compact radio sources may 
also be correlated. In this case the scatterer and the 
receiving antenna act like an adding interferometer to 
produce ripples in the spectrum. However, these effects 
are extremely small since a 1 Jy source results in un- 
der 1 mK of antenna temperature for the dipole-based 
EDGES system and the reflected signal is much smaller 
still. The ground reflection has been eliminated by plac- 
ing the antenna on the ground. A brief discussion on the 
impact of these effec ts in rad io astrono my m easurements 
can be found also in iRohlfs fc Wilsod ()1996D . 

3.2. Hardware Configuration 

Following a careful choice of the deployment site for 
the experiment, the remaining sources of systematic un- 
certainty result from the hardware design of the system. 
The EDGES system consists of three primary modules: 
1) an antenna, 2) an amplifier and comparison switch- 
ing module, and 3) an analog-to-digital conversion and 
storage unit. The antenna, shown in Figure [l] is a "fat" 
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di pole-based design derive d from the four-point antenna 
of ISuh et all (|2003l 120041 ). The design was chosen for 
its simpHcity and its relatively broad frequency response 
that spans approximately an octave. The response of the 
antenna was tuned to 100 to 200 MHz by careful selection 
of the dipole dimensions. In order to eliminate reflections 
from the ground and to reduce gain toward the hori- 
zon, the antenna is placed over a conducting mesh that 
rests directly on the ground. The mesh is constructed 
from thin, perforated metal sheets to reduce weight and 
is shaped to match an octagonal support structure below 
the ground screen. The diameter of this ground screen is 
approximately 2 m. 

Although the antenna is constructed with perpendicu- 
lar dipoles capable of receiving dual linear polarizations, 
only one polarization of the crossed-dipole is sampled 
by the receiver in order to reduce the cost of the sys- 
tem. This is acceptable since the spatially averaged all- 
sky spectrum is expected to have essentially no polarized 
component. The Galactic foreground does exhibit strong 
polarization in certain regions, such as the "fan region" 
around ^ « 140°, 6 « 8°, which has an extend ed polar- 
ized component of about 3 K (( WilkinsonI 1 1 9731 ) . Such a 
region could produce a ripple in the measured spectrum 
from a single linear polarization as the polarization angle 
rotates with frequency. Under the worst circumstances, 
if such a region were located at the peak of the EDGES 
beam, the magnitude of the ripple could reach ~ 50 mK. 
Away from the Galactic plane, however, where EDGES 
observations are generally targeted in order to reduce 
the system temperature, it is more likely that the effects 
of polarization would be at least an order of magnitude 
lower. Furthermore, if the rotation measure (RM) is of 
order 10 rad m~^, then the polarized component could be 
averaged out over ~ 1 MHz. Nevertheless, in future ver- 
sions of EDGES, both ports of the antenna will be sam- 
pled in order to check for polarization effects and other 
systematic effects that result from the non-uniformity of 
Galactic radiation. 

A dipole antenna is naturally a balanced electrical sys- 
tem. To convert from the balanced antenna leads to 
the unbalanced receiver system (in which one lead is 
grounded), a short coaxial cable enclosed in a clamp- 
on split ferrite core with a high impedance is used as a 
common-mode choke balun^ and is connected directly to 
the terminals of the antenna with the central conductor 
fastened to one element and the braided shielding to the 
other. 

The amplifier module consists of two stages that are 
contained in separate aluminum enclosures to reduce cou- 
pling between the low-noise amplifiers. Each stage pro- 
vides 33 dB of gain for a total of 66 dB. Bandpass filtering 
of the signal is also performed in the second stage, and 
the resulting half-power bandwidth spans approximately 
50 to 330 MHz. The amplifier chain can be connected 
through a voltage controlled three-position switch to one 
of three inputs: the antenna, an ambient load, or an 
ambient load plus a calibration noise source. Switching 

* This balun provides a 1:1 impedance transition and operates on 
th e same princ iple as the quarter wavelength sleeve balun described 
bv lKraua Ill988l . page 742). The ferrite provides a high impedance 
over a wide frequency range to reduce the common-mode currents, 
whereas the sleeve balun provides a high impedance over only a 
limited frequency range close to the quarter wavelength resonance. 



between the ambient load and the antenna provides a 
comparison to subtract spurious instrumental signals in 
the measured sky spectrum. 

Impedance mismatch between the antenna and the am- 
plifiers causes reflections of the sky noise within the elec- 
trical path of the instrument that produce an undesir- 
able sinusoidal ripple in the measured spectrum due to 
the frequency-dependence of the phase of the reflections 
at the input to the amplifier. To reduce the effects of 
these reflections in EOR measurements, the input to the 
amplifier chain is connected directly to the balun on the 
antenna (with no intermediate transmission cable), as 
shown in Figure [1] While absolute calibration is limited 
in this configuration by the effect of the unknown phases 
of the reflections on the measured spectrum, the com- 
pact size of the antenna and the small signal path delays 
result in a smooth spectral response. 

The amplifier module is connected to the analog-to- 
digital conversion module by three low-loss coaxial trans- 
mission cables. The cables provide power, switching 
control, and signal transmission, respectively. Common- 
mode current on these cables (i.e. current that is on the 
outer surface of the shielding in the coaxial cable, or cur- 
rent that is unidirectional on both the central conductor 
and inner surface of the shielding) is also capable of pro- 
ducing refiections and additional sinusoidal ripples in the 
measured spectrum. The ferrite core balun used between 
the antenna and amplifiers allows common-mode current 
of approximately 10% of the differential mode. Although 
most of this current is transferred to the ground screen 
by direct contact between the amplifier module casing 
and the ground screen, some current persists and leaks 
through the casing of the amplifier module and onto the 
shielding of the three cables connecting the amplification 
module to the analog-to-digital conversion module. To 
reduce this current to less than 0.005% of the differen- 
tial mode current, additional clamp-on ferrite cores are 
placed every meter on the transmission cables. 

Finally, the analog-to-digital conversion is accom- 
plished with an Acqiris AC240''^ 8-bit digitizer with max- 
imum dynamic range of 48 dB (although, in practice, 
the effective dynamic range was substantially lower due 
to coupling between the digital output of the converter 
and its input). The AC240 uses an embedded field pro- 
grammable gate array (FPGA) to perform onboard Fast 
Fourier Transform (EFT) and integration of the power 
spectrum in realtime. The spectrometer is clocked at 
1 GS/s and the Fourier transform processes 16,384 chan- 
nels, giving a bandwidth of 500 MHz and a raw spectral 
resolution of about 30 kHz. T he broadband spec trome- 
ter employs the FPGA code of iBenz et al.l (|2005f l and a 
Blackman-Harris window function is used to improve the 
isolation between neighboring frequency channels at the 
expense of reducing the effective spectral resolution to 
122 kHz. The unit is contained on a CompactPCI card 
connected to a host computer. The digitizer and host 
computer, along with a power transformer and interface 
circuitry for controlling the amplifier module with the 
serial port of the computer, are enclosed in an aluminum 
box to prevent self interference. 

^ http:/ /www. acqiris. com/products/analyzers/cpci-signal- 

analyzers/ac240-platform.html 
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Fig. 3. — Integrated spectrum used for upper limit analysis 
of reionization signal. The sky temperature, Tgi^y, is an estimate 
based on modeled values for cable losses and no correction for an- 
tenna reflections. The spectrum represents the best 10% of the 
data from observations over two nights. It is selected by discard- 
ing individual observation cycles (see § 13.31 1 containing periods of 
particularly intense radio frequency interference. A total of ap- 
proximately 1.5 h of integration is included (3.75 h including the 
ambient load and calibrator noise source measurements in each cy- 
cle). The black curve shows the spectrum after de- weighting the 
interferers (shown in gray) present in the retained observations. 

3.3. The Measured Spectrum 

The measured spectra from each of the three switch 
positions can be combined to produce a cahbrated esti- 
mate of the true sky spectrum. The three spectra are 
given by 

po = 9 (Tl+Tr) (1+no) 

Pi = g {Tl + Tr + Teal) (1 + ni) (2) 

P2 = g (Ta + Tr) (1 + 7^2) 

where the explicit frequency dependence of each term 
has been dropped and pq is the spectrum for the am- 
bient load, pi is the spectrum for the ambient load plus 
calibration noise, and p2 is the spectrum for the antenna. 
In this terminology, g is the gain, Tl is the ambient load 
temperature, Tr is the receiver noise temperature, Tcai is 
the calibration noise temperature, and Ta is the antenna 
temperature. Thermal uncertainty in the measurements 
is explicitly included in the Gaussian random variables 
riQ, ni, and 712, the magnitudes of which are given by 
rii — [e b Ti)~^l'^ ^ where e — 0.5 is the efficiency for 
the Blackman-Harris window function (which could be 
improved to 0.93 by processing two sets of overlapping 
windows), h — 122 x 10^ Hz is the resolution bandwidth, 
and Ti is the integration time in seconds (for each switch 
position, i). Temporarily setting the noise terms to zero, 
rii — > 0, and solving for the antenna temperature yields 

TA=T,al^^^+TL. (3) 

Pi -Po 

In practice, the impedance match between the antenna 
and receiver is not perfect and some of the incident sky 
noise may be reflected back out of the system. This 
produces deviations between the derived sky tempera- 
ture found using Equation [3] and the true sky spectrum. 
Independent measurements of the impedance mismatch 
can be used to correct these deviations by applying a 
frequency-dependent multiplicative factor to Ta that is 



proportional to the inverse of the reflection coefficient. 
For the EDGES system, this correction was measured 
by two methods: we used a network analyzer in the lab- 
oratory to determine the impedance of the antenna, and 
we reconfigured the system in the field with a long cable 
inserted between the antenna and amplifier module so 
that reflections between the two elements were visible in 
the measured spectrum and could be used to calibrate 
the reflection coefficient. In both sets of measurements, 
the corrections were found to be small (of order 1%) and 
smooth (able to be fit by a low-order polynomial in fre- 
quency) over the band of interest. For the remainder 
of this paper, we will ignore this correction since its ef- 
fects are easily absorbed by the polynomial fit technique 
used to constrain the redshifted 21 cm contribution to 
the spectrum. 

Adding the noise terms back in and solving in the limit 
that Teal ^ {Tl « Ta) > Tr results in an estimate of 
the thermal uncertainty per frequency channel of approx- 
imately 

ATA,rms ~ \JnliTL + Tr^ + uHTl)^ -f nl{TA + Tr^. 

(4) 

For optimal efficiency, the three terms contributing to 
ATa rms should bc comparable in magnitude. Substi- 
tuting Tl = 300 K, Ta = 250 K and T^^ = 20 K, we find 
that the terms are comparable as long as approximately 
equal time is spent in each switch position. In addition, 
a 1 hour integration in each switch position (3 hours to- 
tal) will result in a thermal uncertainty in the estimate 
of the antenna temperature of AT/i,rms ~ 35 mK within 
each 122 kHz frequency channel. 

To acquire a series of estimates of the true sky spec- 
trum using this technique, software on the host computer 
cycles the amplifier module between the three switch po- 
sitions and triggers the digitizer to acquire, Fourier trans- 
form, and accumulate data for a predefined duration at 
each of the switch positions. The integration durations 
per switch position are T{o.i,2} = {10,5,10} seconds for 
the ambient load, ambient load plus calibration noise 
source, and antenna, respectively, giving a duty cycle 
of about 40% on the antenna. This loop is repeated ap- 
proximately every 25 seconds for the duration of the ob- 
servations and the resulting measurements are recorded 
to disk. 

4. INITIAL RESULTS 

The EDGES system was deployed at the radio-quiet 
Mileura Station in Western Australia from 29 Novem- 
ber through 8 December 2006. These dates were chosen 
such that the Galactic center would be below the hori- 
zon during most of the night, keeping the system tem- 
perature as low as possible for the measurements. The 
system was located approximately 100 m from the near- 
est buildings in a clearing with no nearby objects and no 
obstructions above ~ 5° on the horizon, and the antenna 
was aligned in an approximately north-south/east- west 
configuration. The system was operated on 8 consecu- 
tive nights during the deployment, with 5 of the nights 
dedicated to EOR observing. In total, over 30 h of rel- 
evant drift scans were obtained, but strong, intermit- 
tent interference from satellites complicated the mea- 
surements and only approximately 8 h of high-quality 
observations were retained as the primary data set. Al- 
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Fig. 4. — Residuals after subtraction of seventh-order polynomial 
fit to measured spectrum shown in Figure |3] The gray line is the 
raw spectrum with 122 kHz resolution. The black line is after 
smoothing to 2.5 MHz resolution to reduce the thermal noise to 
below the systematic noise. The rms of the smoothed fluctuations 
is approximately 75 mK (see Figure |5)l. 

though the sateUite interferers that complicated the mea- 
surements were narrow-band and, in many cases, were 
easily removed through excision of the effected spectral 
channels, the limited dynamic range of the EDGES sys- 
tem resulted in clipping of the analog-to-digital converter 
and corruption of the full band during especially strong 
transmissions. This required all channels to be omitted 
from the data set in those instances. In particular, it 
was found that the low Earth orbiting satellites of Orb- 
comm (transmitting between approximately 136 and 138 
MHz), as well as satellite beacons (at 150 MHz) from 
discarded spacecraft were particularly troublesome. The 
Orbcomm activity was somewhat variable and usually 
decreased during the night. The typical duration of a 
pass was approximately 15 minutes, during which time 
the power in the satellite signal could easily reach an or- 
der of magnitude greater than the integrated sky noise 
over the band. While previous ob servations at the site 
with prototype MWA equipment (|Bowman et al.l [20071 ) 
have demonstrated (in a subset of the full target band) 
that it is possible to reach the sensitivities required for 
EDGES despite the satellites and other sources of inter- 
ference, improvements to the EDGES digitizing system, 
such as an upgrade of the analog-to-digital converter to 
10 or 12 bits, would certainly help to alleviate the diffi- 
culties encountered during this observing campaign and 
increase the usable fraction of measurements. 

From the primary data set remaining after transient 
RFI exclusion, a stringent filter was applied to select the 
best 1.5 h of sky-time when transient interference signals 
were weakest. The final cut of data included measure- 
ments from multiple nights and spanned a range of local 
apparent sidereal time (LST) between and 5 h. The sky 
temperature at 150 MHz derived from the system during 
this period was found to have a minimum of ~ 240 K at 
about 3 h LST and a maximum of ~ 280 K at 5 h LST. 
The integrated spectrum generated from these measure- 
ments is shown in Figure [3l Frequency channels contain- 
ing RFI were identified in the integrated spectrum by an 
algorithm that employs a sliding local second-order poly- 
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Fig. 5. — Characteristic amplitude of the residuals to the poly- 
nomial fit as a function of integration time on the sky. The rms 
follows a thermal {br)"^/'^ dependency until saturating at a con- 
stant 75 mK noise level due to the instrumental errors introduced 
into the measured spectrum. The dotted lines are guides for the 
eye showing a (fer)"^/^ profile and a constant 75 mK contribution. 

nomial fit and iteratively removes channels with large er- 
rors until the fit converges. The affected channels were 
then weighted to zero in subsequent analysis steps. To 
look for small deviations from the smooth foreground 
spectrum, a seventh-order polynomial was fit to the mea- 
sured spectrum between 130 and 190 MHz (where the 
impedance match between the antenna and receiver was 
nearly ideal) and subtracted. 

The residual deviations in the measured sky spectrum 
after the polynomial fit and subtraction are shown in Fig- 
ure [H The rms level of systematic contributions to the 
measured spectrum was found to be AT^ms ~ 75 mK, a 
factor of 3 larger than the maximum expected red- 
shifted 21 cm feature that would result from a rapid 
reionization. Although it is not obvious in Figure [H the 
variations in the residuals are due to instrumental con- 
tributions and not thermal noise. The large variations 
between 163 and 170 MHz are due to the 166 MHz PCI- 
bus clock of the AC240 and computer, while the gap cen- 
tered at approximately 137 MHz is due to RFI excision of 
the Orbcomm satellite transmissions over a region span- 
ning more than 2.5 MHz. Analysis of the dependence of 
ATrms on integration duration is shown in Figure [5] and 
illustrates that the rms of the residuals follows a ther- 
mal profile ~ (5r)~^/^ initially and then saturates to a 
constant value. After smoothing to 2.5 MHz resolution 
(Az « 0.2), the instrumentally dominated 75 mK thresh- 
old is reached in approximately 20 minutes (1200 s) of 
integration on the sky (50 minutes of total integration 
in all three switch positions). Reordering the individual 
25-second observation cycles used in the full integration 
does not change the behavior in Figure [5l and longer in- 
tegrations (up to approximately 3 h of sky time) , using 
observation cycles with more intense interference, contin- 
ued to decrease the thermal noise, but leave the spurious 
signals and systematic effects unchanged. 

4.1. Limits on Reionization History 

Although the sensitivity level of the initial observations 
with the EDGES system was limited by instrumental ef- 
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fects in the measured spectrum at a level greater than 
the expected maximum contribution due to redshifted 
21 cm emission, weak constraints can still be placed on 
the reionization history. In addition, it is possible to 
make a quantitative assessment of how much improve- 
ment must be made before significant constraints are 
possible, as well as to characterize the best-case outcome 
of future efforts using similar approaches. To begin, we 
introduce a model for the sky spectrum such that 



Tskyiv) = Tgal[v) + T„nb + T2l{v) 



(5) 



where Tgai represents the contribution of all the fore- 
grounds (and is dominated by the Galactic synchrotron 
radiation), Tcmb = 2.73 K is the CMB contribution, and 
T21 is the specific form for the frequency-dependence of 
the redshifted 21 cm emission during the transition from 
the fully neutral to fully ionized IGM. This model ne- 
glects any directional or temporal variation in T^ky and, 
therefore, implicitly assumes an angular average over the 
antenna beam and a time average over the drift scan mea- 
surements performed for the experiment. Since Tcmb and 
T21 are taken to be constant over the sky, only the Tgai 
contribution is affected by this simplification. This does 
not impact the result, however, as long as the foreground 
emission varies slowly on the sky and the antenna pat- 
tern changes slowly with frequency — conditions that are 
presumed to be met in the high Galactic latitude region 
sampled by the dipole-based EDGES system. As a test 
of this assumption, we calculated the residuals after the 
polynomial fit for a bright source with flux comparable to 
Gas A (1400 Jy at 100 MHz) and spectral index (3 = 2.77 
at various positions in the antenna beam using simulated 
beam patterns to determine the frequency-dependence. 
We found, in all cases, less than a ~ 50 /iK residual. 

During the reionization epoch, we define T21 to be 
given by 



r2i(z) = Ar2ii<i 1 + cos 



7r(Zr 



Az/2) 



Az 



(6) 



where AT21 is constant and is the maximum amplitude of 
the redshifted 21 cm contribution, is the redshift when 
XHi{zr) = 0.5, Az is the total duration of the reioniza- 
tion epoch, and we use v — 1420/ (1 -f z) MHz to convert 
back to frequency units. Before the reionization epoch 
[z > Zr -\- Az/2), T21 = AT21, while after reionization 
[z < Zr — Az/2), T21 = 0. The exact form of the tran- 
sition used for T21 has little influence on the outcome of 
the constraints as long as it is reasonably smooth. Fig- 
ure [2] illustrates the modelled redshifted 21 cm spectrum. 
The free parameters in the model are z^, Az, and AT'21. 

For the EDGES best-response frequency range, a cen- 
ter redshift around z^ = 8 allows the largest range of 
Az to be explored. By simulating the combined sky 
spectrum, Tgky, for a range of the two remaining free 
parameters, we can determine the rms of the residuals 
that would remain following the polynomial fit used in 
the EDGES data analysis. Gomparing the rms of the 
residuals in the models to the 75 mK rms of the initial 
measurements gives a good estimate of the region of pa- 
rameter space ruled out so far. Figure [5] illustrates the 
results of this process. The line defining the ruled-out re- 
gion is computed by finding the locus of parameters that 
make the rms residuals in the model equal to 75 mK. 
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Fig. 6. — Constraints placed by EDGES on the redshifted 21 cm 
contribution to the sky spectrum. The dark-gray region at the 
top-left is the portion of the parameter space ruled out by the 
initial EDGES results with ATrms = 75 mK (solid line). The 
dashed lined labelled ATrms = 7.5 mK and the dotted line labelled 
ATrms = 3 mK indicated the constraints that could be placed on 
reionization if the experimental systematics were lowered to the 
respective values. The light-gray region along the bottom is the 
general range of parameters believed to be viable. The redshifted 
21 cm contribu tion to the spectrum is modelled according to the 
description in ii l4.1l with Zr = 8. 

While a more statistically robust analysis is clearly pos- 
sible, little benefit would be gained for the initial mea- 
surements presented here due to the severe systematic 
effects present in the spectrum. 

From Figure [6l it is clear that the initial results con- 
strain only a small portion of parameter space that is well 
outside the expected region for both the intensity of the 
redshifted 21 cm signal and the duration of reionization. 
The best constraint, in the case of a nearly instantaneous 
reionization, is that the redshifted 21 cm contribution to 
the spectrum is not greater than about AT21 < 450 mK 
before the transition. Reducing the systematic contribu- 
tions in the measured spectrum by more than an order 
of magnitude to Ar^ms < 7.5 mK would begin to allow 
meaningful constraints, while an improvement of a fac- 
tor of 25 to ATj-ms « 3 mK would be able to rule out 
a significant portion of the viable parameter range and 
constrain Az > 2. In principle, such an improvement is 
possible with minor modifications to the EDGES system. 
Reaching a systematic uncertainty below ^ 3 mK, how- 
ever, is likely to be infeasible without a redesign of the 
experimental approach because errors in the polynomial 
fit to the overall power-law-likc shape of the sky spec- 
trum, Tskyiy), are the dominant source of uncertainty 
below that level in the current approach. The sharp 
cut-off at Az « 2 in parameter space is the result of 
using a seventh-order polynomial to fit a 60 MHz band- 
width, thus yielding a maximum residual scale size of 
order 10 MHz, which corresponds to Az « 2 at z « 8. If 
the same polynomial could be reasonably fit to a larger 
bandwidth (or a lower-order polynomial fit to the ex- 
isting bandwidth), then Az could be probed to larger 
values. 

5. CONCLUSION 



9 



In principle, useful measurements of the redshiftcd 
21 cm background can be carried out with a small radio 
telescope. These measurements would be fundamental to 
understanding the evolution of the IGM and the EOR. In 
particular, the global evolution of the mean spin temper- 
ature and mean ionization fraction of neutral hydrogen 
in the high redshift IGM could be constrained by very 
compact instruments employing individual radio anten- 
nas. We have reported preliminary results to probe the 
reionization epoch based on this approach from the first 
observing campaign with the EDGES system. These ob- 
servations were limited by systematic effects that were 
an order of magnitude larger than the anticipated sig- 
nal and, thus, ruled out only an already unlikely range 
of parameter space for the differential amplitude of the 
redshifted 21 cm brightness temperature and for the du- 
ration of reionization. Nevertheless, the results of this 
experiment indicate the viability of the simple global 
spectrum approach. 

Building on the experiences of these initial efforts, 
modifications to the EDGES system are underway to 
reduce the residual systematic contribution in the mea- 
sured spectrum and to expand the frequency coverage of 



the system down to 50 MHz or lower in order to place 
constraints on the anticipated transition of the hyperfine 
line from absorption to emission as the IGM warms be- 
fore the EOR. Constraining the redshift and intensity of 
this feature would be very valuable for understanding the 
heating history of the IGM and, since the transition has 
the potential to produce a step-like feature in the red- 
shifted 21 cm spectrum with a magnitude over 100 mK 
(up to a factor of 4 larger than the amplitude of the step 
during the reionization epoch), it may be easier to iden- 
tify than the transition from reionization — although the 
sky noise temperature due to the Galactic synchrotron 
foreground increases significantly at the lower frequen- 
cies, as well. Through these and other global spectrum 
efforts, the first contribution to cosmic reionization sci- 
ence from measurements of the redshifted 21 cm back- 
ground will hopefully be achieved in the near future. 
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